Electrodynamics of Amorphous Media at Low Temperatures 
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Amorphous solids exhibit intrinsic, local structural transitions, that give rise to the well known 
quantum-mechanical two-level systems at low temperatures. We explain the microscopic origin of 
the electric dipole moment of these two-level systems: The dipole emerges as a result of polarization 
fluctuations between near degenerate local configurations, which have nearly frozen in at the glass 
transition. An estimate of the dipole's magnitude, based on the random first order transition theory, 
is obtained and is found to be consistent with experiment. The interaction between the dipoles is 
estimated and is shown to contribute significantly to the Griineisen parameter anomaly in low T 
glasses. In completely amorphous media, the dipole moments are expected to be modest in size 
despite their collective origin. In partially crystalline materials, however, very large dipoles may 
arise, possibly explaining the findings of Bauer and Kador, J. Chem. Phys. 118, 9069 (2003). 
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I. INTRODUCTION 

Glasses are frozen liquids and thus lack long-range or- 
der, yet the differences in material properties between 
amorphous materials and crystals are often rather sub- 
tle. Crystalline samples themselves are rarely flawless 
and thus contain a number of imperfections such as point 
defects, dislocations or grain boundaries of various sorts. 
These tend to further mask the difference. The size of de- 
fects in crystals ranges over many scales, while in glasses, 
the static heterogeneity in the atomic arrangement ap- 
pears comparable to the molecular size itself. Simple 
molecular glasses thus seem perfect candidates for de- 
scription as isotropic continuum, at long enough wave 
lengths. For instance at cryogenic temperatures, when 
the de Broiglic wave-length of a thermal phonon at ~ IK 
exceeds the lattice spacing by three orders of magnitude 
or so, continuum theory would be thought to hold to high 
accuracy. Yet surprisingly, there clearly exist degrees of 
freedom numbering in great excess of the Debye density 
of states, leading to extra heat capacity and phonon scat- 
tering in all amorphous materials Q. Here we examine 
the electrodynamics of these degrees of freedom. 

Since Rayleigh scattering is too weak to account for 
the observed magnitude of sound attenuation in glasses, 
internal resonances must be involved, in the form of (in- 
harmonic structural rearrangments, in order to explain 
the data. The well known, empirical two-level system 
(TLS) theory presumes such resonances exist 0, 01 ■ Sim- 
ply postulating a flat energy spectrum and a frequency 
independent coupling to the phonons accounts for all 
the gross features of the low T anomalies (for reviews, 
see 0, |M @)- Direct microscopic evidence of the two- 
level nature of such entities comes both from the phonon 
echo experiments Q and, relatively recently, from single- 
molecule experiments at cryogenic temperatures (see e.g. 



At these temperatures, the TLS picture is internally 
consistent in so far as the structural transitions (ST) can 
be defined as local, and thus, tautologically, sufficiently 
weakly interacting. One may therefore speak of a mul- 
tilevel system at the location of each transition whose 
behavior reduces to a TLS behavior at low enough T. 
We may call this a tunneling center (TC). In 1986, Free- 
man and Anderson |9j showed that the magnitude of the 
TLS density of states is apparently correlated with the 
phonon coupling. This results in a universality of the ra- 
tio of the phonon mean free path Z m f p to its wave-length 
A: Zmfp/A ~ 150, for all insulating glasses at T < IK. 
This universality seems hardly coincidental |10| , however 
understanding the origin of the universality requires a mi- 
croscopic picture of molecular motions in glasses. (The 
large size of the factor ~150, too, was a puzzle [l0|.) 

The Random First Order Transition (RFOT) Theory 
of the glass transition 0, 0, 0, 0] provides an 
appropriate microscopic picture of the motions in glass. 
Most commercial and laboratory glasses are made by 
quenching supercooled melts. In the deeply supercooled 
regime, most liquid motions are activated transitions be- 
tween distinct aperiodic states of comparable energy, dur- 
ing which the current structural arrangement in a local 
region is replaced by another, quite different arrangement 
that nevertheless fits its environment. The size of the 
reconfigurable region, £, grows with decreasing temper- 
ature, and reaches about 5-6 molecular units across by 
T g , i.e. the glass transition temperature corresponding to 
the 1 hour time scale. This size is predicted to be univer- 
sal, within logarithmic accuracy, for all substances. At 
any point in time, above T g , the liquid can be thought 
of as a mosaic of such cooperative regions ^4|, most of 
nearly the same size, but otherwise with distributed bar- 
rier heights and transition energies. Upon freezing, a par- 
ticular mosaic pattern sets in and undergoes relatively 
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slow changes, called aging. The aging speed depends 
on the quench depth [16j. A sufficient fraction of the 
structural transitions have small enough energy change 
and barriers so that when they occur, they can account 
for the cryogenic anomalies, some of which were men- 
tioned above: the density of states and the universality 
of phonon scattering > the Boson Peak , but also 
the anomalous Griineisen parameter, the so called "fast" 
TLS systems and more HlUl. According to the RFOT 
theory, the universality of the l m f p / A ratio directly follows 
from the universal cooperative region size (£/a) 3 ~ 200 
at the glass transition temperature T g , where a is the 
molecular length scale. During a structural transition, a 
relatively large, ~ 200, compact set of small units moves 
in a stage-wise fashion. This corresponds to the motion 
of the domain wall, which separates the two alternative 
arrangements, through the compact region. At cryogenic 
temperatures, these motions occur by tunneling. Despite 
their collective nature, such tunneling events are possible 
because of the enormous multiplicity of alternative struc- 
tural states and of low-barrier paths connecting pairs of 
states: an amorphous sample actually resides in a high 
energy density state, well above its lowest energy, perfect 
crystalline state. Consistent with the facility of tunnel- 
ing is the smallness of individual atomic displacements 
during each transition. Their amplitude is roughly equal 
to the Lindemann length d^. This length is typically one 
tenth of the characteristic lattice spacing a and is nearly 
the same for all substances. The precise identity of the 
"molecular unit" , or "bead" depends on the specific sub- 
stance, but usually corresponds to a few atoms. (See 
for a detailed discussion.) 

In this article, we use the microscopic picture of the 
two-level systems provided by the RFOT theory to es- 
timate the coupling of the transitions to external elec- 
tric fields. Clearly, such a coupling must be present be- 
cause individual molecular bonds, that possess electric 
dipoles, rotate during transitions. Since these couplings 
directly enter into spectral hole-burning experiments j2J 
and can also be directly probed in single molecule exper- 
iments 0, it is important to know how much collective 
excitations interact electrodynamically with probes and 
external electric fields. 



II. INTERACTION OF A SINGLE TUNNELING 
CENTER WITH EXTERNAL FIELD 

A. Many-Body Origin of the Transition-Induced 
Dipole Moment 

To set the stage, let us briefly review the assumptions 
of the traditional molecular models of dielectric response 
of insulating media. One often assigns an electric dipole 
value to an individual molecule, or to a molecular bond 
connecting distinct atoms in a condensed phase. In a 
dilute liquid made up of polar molecules, the medium 
polarizes in a field since the dipoles prefer to orient along 



the field's direction at the cost of losing their rotational 
freedom. Even without permanent electric dipoles, a di- 
electric response occurs due to polarizability: An exter- 
nal field mixes in higher energy molecular orbitals, which 
generally lack inversion symmetry. Classically, this quan- 
tum mechanical response can be imitated as two harmon- 
ically bound opposite charges that separate after a field 
is turned on. In a polar substance, this polarizability 
also changes the length or orientation of the permanent 
dipoles. 

While in a fluid the dipoles can freely reorient, dur- 
ing a structural transition (ST) in glass, the dipole is re- 
strained: Each individual bead within a tunneling center, 
or "domain" , moves only about the Lindemann length 
d,L, as illustrated in Fig^ Suppose, for the sake of ar- 
gument, one can break up the set of all the beads within 
the domain into distinct pairs. During a transition, each 
such pair - and hence the corresponding "bond" - rotates 
about cIl/cl ~ 0.1 radian. The amorphous lattice gener- 
ally exhibits no symmetry. There is therefore, typically, 
excess charge, however small, on each atom. Assume 
the effective individual charges remain the same during 
such a transition. One can thus unambiguously assign a 
permanent, point-like electric dipole to each "bond" in- 
troduced above. As schematically shown in Fig^ a total 
dipole moment, /xt = Yli ^Mi: may be generated during 
a structural transition, that would couple to an external 
electric field E with energy —^E. 
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FIG. 1: Shown on the left is a fragment of the mosaic of 
cooperatively reconfiguring regions in the supercooled liquid, 
with a denoting the lattice spacing (more precisely "bead" 
spacing). £ is the cooperative region size; cLl is the typical 
bead displacement during a transition. (The shown magni- 
tude of £ corresponds to a temperature near T g on 1 hour 
scale.) The two sets of circles - solid and dashed ones - de- 
note two alternative structural states. The expanded portion 
shows how rotation of a bond leads to generating an elemental 
dipole change during a transition, where the partial charges 
on the two beads are ±£g. 

In molecular glasses, the bond dipoles are fairly easy 
to assign. Generally, the assignment of point-like dipole 
moments to individual bonds is, strictly speaking, non- 
unique. It would be rather difficult in the case of a highly 
networked, covalently bonded substance, such as amor- 
phous silica, but as we shall see, other arguments for 
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such systems give similar results. These arguments use 
the measurable piezoelectric properties of corresponding 
crystals to unambiguously extract the coupling to fields. 
In weakly bonded molecular glasses held together by Van 
der Waals forces, the point like dipole view is already a 
good approximation. 

With this in mind, an order of magnitude estimate of 
the dipole moment of a tunneling center can be made: 
The Coulomb charge on a bead does not exceed a frac- 
tion £ < 1 of an elementary charge q, which is close in 
magnitude to the electron charge e: q ~ e. An individ- 
ual electric dipole change is therefore A/j, ~ CMmoif^i/a), 
where {d^/a) is the rotation angle, as already discussed, 
and /imoi = (qa is the elemental dipole magnitude associ- 
ated with each bond. The number of pairs that reorient 
in a structural transition is iVd = (£/ a ) 3 /2- The (At 
distribution is, of course, centered at the origin. Since 
the dipole forces are a small part of the energetics of the 
glass transition, the dipole motions are expected to be 
only weakly connected with each other. Therefore the 
individual dipoles A/2j make up a "random walk" of 
steps (in 3D). As a result, the generic value of dipole 
change for the transition is given by the width of the 
total displacements during such a walk: 

Mr ~C(go)[«/a) 3 /2] 1/2 (d £ /a), (1) 

If the elementary dipole rotations are correlated, one may 
introduce an additional factor - like the Kirkwood g fac- 
tor of liquid theory g is typically of order 2. At 
first we might imagine a great deal of variability for the 
quantity [(£/ a) 3 /2] 1 / 2 (d^ / a) but in fact it is nearly uni- 
versally (!) equal to unity for all substances, a is typi- 
cally a couple of angstroms, implying qa corresponds to 
~ 10 Debye. (eA <-> 4.8 Debye.) £ is expected to be 
well less than unity, with 0.1 or less being a reasonable 
generic estimate. We thus obtain that fix is of the order 
1 Debye or less, consistent with experiment. Note the 
magnitude of fiT is rather modest - only of order the size 
of a typical individual dipole moment - despite the large 
number of particles constituting a domain. The relative 
smallncss of the two- level system dipole moment is due to 
the small deflection angle ~ 0.1, and the apparent small- 
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where e cy ±{u>) is the dielectric susceptibility outside the 
domain. The dielectric constant inside, eTLs(^), con- 
tains both the usual (high frequency) polarizability of 
the material and the polarization due to the dipolar dis- 
placements accompanying the transition. Even though 



ness of the partial charge (q. There are deep reasons for 
both: The former stems from the particular magnitude 
of atomic displacements during a transition: it is equal 
to the typical thermal displacement at the mechanical 
stability edge, i.e. the Lindemann length [l4|. The latter 
is probably related to the intrinsic difficulty of making 
ionically bonded aperiodic structures, which imposes an 
upper bound on the value of £; this will be discussed in 
due time. 

The sub-linear scaling of the dipole moment with the 
domain volume £ 3 , in Eq.Q, is worth noting: The tun- 
neling transition dipole moment is not a bulk response, 
and its relation to the material's average bulk dielectric 
properties, as encoded e.g. in the substance's dielectric 
susceptibility e(w) in the fluid phase, is not immediately 
obvious. While it is ultimately the deflections of the 
same elemental dipoles that both give rise to e(u>) and 
the dipole moment, the causes of the deflections differ. 
In contrast with the bulk polarizability, where the de- 
flection magnitude is proportional to the field, the tran- 
sition induced deflections are intrinsic and correspond to 
distinct local structural states. (To illustrate this dis- 
tinction further, we point out that a transition can be 
induced by things other than an AC (!) electromag- 
netic field - a thermal phonon for instance.) Above T g , 
the distinct structural states, that evolve into TLS at 
cryogenic temperatures, are transient metastable struc- 
tures that live typically as long as the a-relaxation time 
14]. The transient structures are (transiently) frozen-in 
elastic fluctuations. Analogously, the intrinsic generated 
dipole moment may be thought of as due to frozen-in 
electric fields - at T g . (A formal connection between gen- 
erated electric field and mechanical stress is discussed in 
the following Subsection.) We can use this line of thought 
to relate the dipole moment to the bulk dielectric prop- 
erties of the material near T g . Suppose the frozen-in field 
is along z direction. The corresponding dipole moment is 
related to this field via an appropriate (frequency depen- 
dent) dielectric constant 6tls(w), but also through the 
fluctuation-dissipation theorem. A self-consistent closure 
of this relationship gives for a spherical region of volume 
£ 3 0: 



(2) 
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the interior and exterior of the domain are chemically 
identical, it is necessary to regard the two e's as distinct, 
since we know a transition occurs within the volume £ 3 . 
A similar relation can be written for the volume occupied 
by a single elemental dipole, with e cx t(w) = eTLs(^) : 



d <Mr„(o)^(t)> dt = [^sn-ij^cxtH + i] ^^ 



kBTg J dt ' ' 47r[2e cxt (o;) + e T Ls(w 
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- t (^A0)^At)) dt 



[e TL sH ~ l][2e TL sM + 1] 
47r[3e T Ls(^)] 



(2a 3 



(3) 
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where the factor (d^/a) 2 on the left reflects that the el- 
emental change in polarization, A/i,, is related to the 
full elemental dipole fii by the rotation angle (d^/a), 
see FigQ] The volume 2a 3 on the right corresponds to 



the volume occupied by a pair of beads, as before. One 
gets, as a result, a frequency dependent generalization of 
Eq.JTJ: 
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,(2), 



M )Wa) 3 /2] 



[2e oxt (w) + l]3e T LsM 



[2e TL s(w) + l][2e ex t(w) + ctlsM] ' 
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where the two-point correlation functions are the t inte- 
grals above. Note ^Mt 3 ( )) = 04) and (vfH )) = 
{ill). Finally, fi t ~ ((qa) = fj, moh as before. 

Note that two adjacent regions are statistically unlikely 
to undergo a structural transition at the same time. The 



physically preferable choice for the external dielectric sus- 
ceptibility e oxt (oj) is therefore the high frequency, elec- 
tronic component of the full dielectric response, which 
we call £oo. With this, Eq.|@J becomes 



A*t ( w ) 



[2eoo + l]3e T Ls(^ 



[2e TL s(w) + l][2ec + eTLs(w)]' 



(5) 



Furthermore, the two-level systems that are active at 
low temperatures correspond to the low barrier side of 
the barrier distribution. This implies one should use the 
eTLs(w - * 00 ) value with regard to the cryogenic phe- 
nomena and, therefore, no extra frequency dependence 
appears in the coupling of the TLS to electric field. As a 
result, we find no significant reaction field correction to 
our earlier argument and obtain 



c.f. Eq.Q. 



B. Piezoelectric View 
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To treat covalent network glasses, where the assign- 
ment of local dipoles or changes is difficult, we now turn 
to a different way to relate the dielectric response, within 
a domain, to frozen mechanical stress. First consider a 
strictly periodic lattice that lacks parity symmetry. Gen- 
erally, the lattice strain induces internal electric fields 
giving rise to piezoelectric behavior. In such a piezoelet- 
ric, the energetics of the strain, in the lowest order, are 



described by the free energy density (see e.g. |24j1 

F = -Kk.lm u ik u lm — T—£ikEiE k + /3j kl^i u kl- (?) 

Here Uik is the standard strain tensor |25|. Xik,im and 
are the stiffness tensor and the dielectric tensor respec- 
tively, and 0i t M is a piezoelectric tensor. (Note, various 
sign conventions and free energies have been used in the 
literature.) The double index summation convention is 
implied throughout, with the exception of letters x, y, 
and z, which will be obvious in the context. Finally, the 
elastic constant A is related to the material's mass density 
p and the speed of sound c s : 



A - pc 2 s 



(8) 



In the absence of external field: D = —Aiv(dF / dE) = 
0, - the internal electric field is simply proportional to 
the strain itself: 



Ei = Aire- 1 ?, 



(9) 



(Note, the j3 tensor has the dimensions of electric field.) If 
expressed in terms of strain only, the free energy density 
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reads, in the absence of external electric field: 

F = \K ,kiUijU ku (10) 

KjM = ^ij,kl + ^ ne mnPm,ijPn,kl- 0- 1 ) 

The total, apparent stiffness, A', can be decomposed 
thereby into a purely "covalent" and a "coulomb" , i.e. 
electrostatic component. The latter contribution is ordi- 
narily quite small, owing to the smallness of the charges 
induced by lattice distortions. Consider a-quartz, for ex- 
ample. Here, only j3 XlXX — 5.2 • 10 4 esu, and /3 x , V z = 
Px.zy = (1/2)1-2 • 10 4 esu, are non-zero 26]. As a re- 
sult, the "coulomb" contribution to the (xx, xx) compo- 
nent of the apparent stiffness tensor is only about one 
percent of the covalent counterpart: X xx , xx — 8.8 • 10 11 
dyne/cm 2 and (47r/ef )/3 2 xx = 7.4- 10 9 dyne/cm 2 . (Here, 
we used e[ = 4.58 [l^J.) The relative size of the A 
vs. (3 magnitudes can be understood as follows: A re- 
flects the energy (density) of the elastic restoring force. 
It is essentially the second derivative of an individual 
atomic potential. Interatomic bonding, be it coulom- 
bic or covalent in character, is ultimately of electrostatic 
origin. One may therefore associate A with the quan- 
tity ^3 Jv(<? 2 / r )lr=a = (<z/a 2 ) 2 , where q is the effective 
charge giving rise to the bond, and the 1/a 3 factor in 
front provides for energy density. (The total first deriva- 
tive of the (full quantum) potential energy is zero, of 
course.) By virtue of being an electric field, (3 roughly 
corresponds to the quantity q' /a 2 , where q' would be the 
partial charge introduced in the previous subsecton. The 
ratio q' / q corresponds, within the present framework, to 
the earlier introduced quantity £. It follows that in a- 
quartz, the partial atomic charge is indeed about one 
tenth elementary charge, since the quantity 

C 2 = (q'/q) 2 ~ /3 2 /A (12) 

is approximately equal to 1/100 in silica, as we just saw. 

Now suppose for a moment that a relation similar to 
Eq.@ exists between the bead displacements within a 
domain and the internal electric field changes generated 
during a transition. (We stress, in an amorphous sam- 
ple such generated field changes are zero, upon spatial 
average, but here we refer to local fields at a particu- 
lar, generally non-centrosymmetric site.) Since dF — 
—DdE /An + . . ., the free energy change in the presence 
of a (small) external field D cx t during the transition is 
given by J y dV AF = -D cxt J v dV AE/4tt, in the low- 
est order in D cx t- (Here we have volume-integrated over 
the reconfiguring domain that correspons with the two- 
level dynamics; D cx t obviously varies sufficiently slowly 
within the domain for realistic frequencies of light, and 
can be taken out of the integral.) The relation of the 
field .Dext to the external field proper depends, of course, 
on the experiment's geometry. We will use an electric 
field E cx t(r) = e _1 £) oxt (r) where e is the average bulk 
dielectric susceptibility (which, of course, is uniform and 



isotropic in an amorphous material). The coupling to 
this field .Eext is consequently given by: 

Ht = ^ J dV AE(r), (13) 

i.e. the generated internal field difference during the tran- 
sition, integrated over the domain. Eci. (|13J) . among other 
things, demonstrates that one may indeed unambigu- 
ously assign a collection of point-like dipoles to the bead 
set within a TC, namely by virtue of the relation E{r) = 
4ircifi(r)n(r), where n(r) is the (coordinate-dependent) 
dipole concentration and the constant c\ ~ 10° should 
be chosen depending on a specific way to incorporate the 
already mentioned cavity effects. In what follows, we will 
outline the microscopic picture of interaction of a tran- 
sition with elastic strain, which will naturally lead us to 
the formula above and the ability to estimate the electric 
moment via material's piezoelectric properties. 

Since individual displacements di during tunneling 
transitions are only one-tenth of the lattice spacing, one 
can indeed describe the corresponding additional elastic 
energy variations, due to the presence of a phonon, by 
a quadratic form of the type from Eq.Q or l|10|) . Fur- 
ther, here one computes relative displacements, not the 
absolute atomic coordinates which are generally difficult 
to calculate. Define fak as the strain tensor due to a 
(long-wave) lattice distortion of a stable lattice. In ad- 
dition, define dik as the "strain" tensor corresponding to 
the set of the tunneling displacements {di}. The full elas- 
tic energy within the domain, given a particular domain 
boundary configuration (call it fib) can be written as: 

F = 2^ij,kl [(4>tj+d l j)((j)ki+d k i)-d i: jd k i}+n{{d lj }, fl b ), 

(14) 

where the energy functional H({dij}, fib) includes all the 
non-linear, many-body interactions giving rise to the ex- 
istence of the many metastable structural minima within 
the domain. The construction of a library of the states 
corresponding to these minima was described in Ref . [Tfij ] . 
The full multiplicity of the local states reveals itself di- 
rectly in calorimetric measurements above the glass tran- 
sition. Here, we are only concerned with the two lowest 
energy states of a region of the otherwise undisturbed 
lattice. The two correspond to the two lowest minima 
of H({dij}, fib)- The size of the domain is, as we have 
seen, chosen so that one is guaranteed to have at least 
one alternative structural state of nearly the same en- 
ergy, and was found to be o nly slightly larger that the 
cooperative region size at T g [17]. Note, by construction, 
the boundary state fib is independent of the phonon field 
(j)ik- The effect of an external mechanical stress on the 
internal displacements within a local, compact region is 
passed on through the boundary, and so the interaction 
of the region with the stress can be expressed through a 
displacement integral over the region surface [25|. The 
cross Xij.M&ijdki term in Ea. (|10|l gives the amount by 
which the energy of a tunneling transition is modified by 
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the presence of a phonon. This therefore gives the TLS- 
phonon coupling. The latter coupling was estimated in 
this way in Ref.0, |23|- A direct computation of the 
4>ij -field-induced change of the transition energy gives: 
AE(<pij) = A'y ki4>ij J v dV dki, where we have integrated 
over the domain volume and taken advantage of the elas- 
tic strain being tfiij nearly constant throughout the do- 
main (A' in the latter equation is the domain-averaged 
value of the atomic force-constant.) The coefficient at the 
<pij gives a (tensorial) coupling of a transition to strain 
according to: 

^TLS,ph = gij&j&z, (15) 

where a z = ±1 is the usual Pauli matrix, and 

9ij = \Ki,U J y dV d kl . (16) 

The volume integral above indeed reduces to a surface 
integral of the tunneling displacements. Consider for ex- 
ample, the term J v dV d xy = J s (d x dx+d y dy)dz etc. The 
coupling to the longitudinal phonons has the most vivid 
form, since da is the divergence of a vector field, i.e. d. 
One gets gu ~ A J s d dS. The atomic displacements at 
T g are typically near the Lindemann length: d ~ dr., - but 
can also be expressed in terms of the elastic constants, 
since the amount of elastic energy contained in a unit cell 
is determined by the temperature itself: a 3 X(d/a) 2 ~ T g . 
(This is by virtue of the fluctuation-dissipation theorem.) 
Estimating the surface integral 0, |2(j introduces addi- 
tional numerical factors and gives, within a factor of two 
or so: 

g ~ y/T gP c*a?, (17) 

where we have used Eq.(|SJ. The result above is easy to 
rationalize on general grounds: Any atomic motions in 
a dense liquid (at T g and otherwise) are either a vibra- 
tion or an anharmonic motion that is part of a struc- 
tural transition. The two excitations must coexist and 
thus be marginally stable against each other, in order for 
both to be present. Such a marginal stability criterion 
gives {(Tx9ij<t>ij) — (^ij,kl<f>ij<f>kl}, as follows from optimiz- 
ing Ea. (|15fl together with the elastic energy ^Xij^kityij^ki 
with respect to fcj , multiplying by and thermally av- 
eraging (at T = T g ). Owing to (o^-) ~ (\<t>ij\), Eq.(T7|) 
follows. 

In the same way that the elastic fluctuations inter- 
act with the atomic displacements, polarization waves 
and external electric sources will interact with the in- 
ternal electric fields generated in the domain during a 
transition. As we have seen, the mechanical response 
characteristics of a transition arise in response to stress 
fluctuations at T g . Analogously we can say, the electric 
moments of the two-level system arise in response to the 
electric field fluctuations at T g . In full analogy with the 
fluctuation-dissipation theorem context, the TLS dipole 
moment will interact with an external field source just as 



it did with the internal electric fields at the moment of 
freezing. We may thus compute the dipole moment by 
substituting the generated electric field from Eq. J3J into 
Eci. (|13|) . bearing in mind that the local lattice and the 
corresponding tensors are no longer subject to any par- 
ticular symmetry. Still, since the lattice locally resembles 
a crystalline lattice, one may choose coordinates, again 
locally, in such a way that the the (3 tensor is maximally 
close to a crystalline one. (Clearly, the m, n sum in 
Eci. (|ll|) is independent of the coordinate choice.) There- 
fore the latter sum will give a comparable result to that 
for a crystal. Summing over the displacement tensor dij 
(contracted with the (5 tensor) is quite analogous to the 
argument leading to Eq.JIJ. As a result one obtains the 
following qualitative estimate: 

MT ^(^a 3 )[(C/a) 3 /2] 1/2 (dL/a), (18) 

where f3 is the local value of the piezoelectric cou- 
pling. (We remind the reader that glasses are on average 
isotropic and thus can not have bulk piezoelectric prop- 
erties. It is only in the frozen state that parity is locally 
broken.) The simple relation 

Cqa ~ (3a 3 (19) 

establishes the connection between the "molecular 
dipole" view of the previous Subsection and the "piezo- 
electric" analysis in this Subsection, c.f. Eci. (|12fl . Us- 
ing the quartz parameters above, one obtains that [At is, 
again, of the order Debye. Note that formula (|18fl uses 
quantities that can be measured independently for sub- 
stances which have a crystalline counterpart. The bead 
size a can be determined from the fusion entropy [l5| . 
The piezoelectric constants are measurable too. The two 
views - one based on molecular moments, the other on lo- 
cal piezoelectricity - are somewhat distinct but are highly 
overlapping: Since the local hybridization pattern on in- 
dividual atoms is intrinsically asymmetric in amorphous 
lattices, partial atomic charges, however small, are al- 
ways expected to be present in glasses, giving rise to 
both local permanent dipoles and local piezoelectricity 
Mixing in a dipolar species would enhance both effects. 
According to Ref.[28||. the two- level systems 's dipole mag- 
nitude is very correlated - nearly proportional - to the 
OH~ ion concentration, in amorphous silica with OH~ 
impurities. Yet extrapolation to small ion concentrations 
shows TLS in silica exhibit an intrinsic dipole moment, 
as was later confirmed by an electric dipole echo study 

C. Electrodynamics and Electroacoustics: 
Connection with Experiment 

In order to discuss experiments of two-level systems 
in glasses, involving external fields, let us first recapit- 
ulate a few aspects of the traditional phenomenological 
description, along with the microscopic explanation. The 
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Hamiltonian of an isolated two-level system, as usually 
written in the low T glass context, is 

Wtls = 2 ecr * + \^° x + 9ij4>ijO- z , (20) 

where e is the transition energy, A is the tunneling ma- 
trix element. (The phonon part of the full Hamiltonian 
is given in Eq-JEJ.) According to Refs. fTrl H^ . the 
TLS that are thermally active at cryogenic temperatures 
have their splitting e distributed according to a simple 
Boltzmann-like law n(e) = ■^-e~ t l T a . This roughly de- 
fines the density of states of the tunneling transitions. 
The density of states, as seen by calorimetry, is time de- 
pendent, because the tunneling matrix elements A are 
widely distributed. According to the semi-classical anal- 
ysis in Refs.mili, thc distribution is P ( A ) K 1/A 1+C , 
where c <C 1 is a small constant (c oc hiojj/kBTg). This 
distribution is close to, but not precisely the same as 
the inverse distribution oc 1/A postulated by the phe- 
nomenological TLS model. Including quantum effects 
reveals that the low splitting two-level systems, i.e. those 
with e« A, are special in the following sense. In such re- 
gions, the excess strain energy of the glass is concentrated 
in the domain wall itself, while the barrier separating the 
two alternative structural states is not high enough to 
keep the domain in any of the classical structural states 
as defined in terms of the classical atomic coordinates. 
Such TLS, with depinned domain walls, give rise to an ex- 
tra piece in the combined P(e, A) distribution 20]. They 
correspond to the so called "fast" two-level systems, in- 
troduced early on phcnomcnologically 30] in order to 
rationalize certain quantitative shortcomings of the orig- 
inal TLS model. The quantum depinning of the domain 
wall has been also called "quantum mixing" j2jj . Finally, 
we note the Hamiltonian in Eq.l|20[l leads to rich relax- 
ational behavior, due to both interaction with phonons 
and the phonon-mediated interaction with other TLS's. 
This has been discussed in detail previously [3l|, ||, 
The effects of the interaction of the dipolc moment 
with a static electric field are actually quite difficult to ob- 
serve under routine laboratory conditions (see e.g. |32|). 
The upper limit for the field is given by the dielectric 
breakdown value and is generically 10 6 V/cm. For the 
typical dipole moment of 0.5 Debye, this implies an inter- 
action energy of only I0 -3 eV ~ 10 Kelvin. The constant 
field strength normally employed is actually an order of 
magnitude weaker, or less. In a field E cx t, the transition 
energy is modified according to e — > (e—fi^E). Typically, 
\HtE\ < IK. This is clearly inferior to the characteris- 
tic energy scale of the TLS spectrum, namely the glass 
transition temperature T g [jjfl • The effect of a constant 
field thus turns out to be very generic because of the in- 
trinsic flatness of the energy distribution: The angular 
part of /^T-Ecxt is, obviously, uniformly distributed re- 
sulting, again, in a flat distribution of thc field-modified 
transition energy. In order to discern such a small en- 
ergy variation, a resonance technique must be employed. 
Just such an experiment was performed by Maier at el. 



[2l|, who took advantage of the possibility to burn very 
narrow holes - only a few MHz - in the chromophorc's 
inhomogeneous spectrum. These authors turn on the 
field immediately after burning the hole and observe the 
(time-dependent) hole broadening, whose overall magni- 
tude depends quadratically on the field strength. Maier 
at el. report the value of \it — 0.4 D for a PMMA matrix. 

It follows from our theory that light and sound couple 
to the two-level system transitions in a very similar way, 
save the dipole character of the TC-photon interaction 
distinct from the tensorial coupling of the transitions to 
the phonons. Indeed, the temperature dependence of the 
speed of light in vitreous silica, as obtained early on in 
by Schickfus at el. [3j|, nearly coincides with the cor- 
responding ultrasonic data [{|. This strongly suggested, 
at the time, that both the electromagnetic and acoustic 
anomalies had the same origin. That the coincidence is 
not purely circumstantial was shown soon afterwords in a 
number of elegant electro-acoustic experiments: Increas- 
ing the AC electromagnetic field leads to saturation of 
the structural transitions and a decrease in ultrasonic at- 
tenuation 0, [3{| . In addition, exposure to the AC field 
affects the acoustic impedance of a glass [35| . Again, the 
sufficient sensitivity of these experiments is due to the 
interaction with the AC field being resonant. Finally we 
mention yet another venue in investigating the TLS cou- 
pling to electric fields, namely the electric dipolar echo, 
see e.g. HQ- 

The dipole moment magnitudes, reported in all these 
experiments on the respective substances, are all of the 
order I Debye, although more recent measurements seem 
to be converging on a fraction of a Debye. Unfortunately, 
the extracted dipole values do not completely agree be- 
tween different experiments. So, for example, Kharlamov 
at el. give relatively low values of 0.2 D and 0.1 D 
for PMMA and PS respectively, based on their field in- 
duced spectral diffusion data. The degree of the quan- 
titative discrepancy is, of course, subject to the detailed 
assumptions on the distribution of the individual TLS 
parameters, various angular averagings etc. 



III. DIPOLE-DIPOLE INTERACTION 

The idea of local structural transitions is internally 
consistent in that the transitions are indeed distinct, 
weakly interacting entities. This is easy to understand 
by considering the moment of vitrification, when a par- 
ticular pattern of mobile regions sets in: A transition 
will be found locally, upon freezing, if at T g it was of 
marginal stability with respect to external mechanical 
perturbation, as delivered by stress waves to the given 
local region. It does not matter, of course, whether the 
source of these waves is thermal elastic fluctuations or 
the other structural transitions. Now, upon having es- 
timated the energy spectrum of the TLS and their cou- 
pling to the phonons |17| . one may check the magnitude 
of the resultant TLS-TLS interaction, mediated by the 
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acoustic waves in the frozen lattice. Such interaction 
self-consistently turns out to be small 20] . We are aware 
of several observable consequence of the interaction. For 
one thing, this interaction tends to quench the sponta- 
neous echo generation 7, 30] , by virtue of dephasing each 
TLS's motion when the TLS precesses about its local 
field. Another, remarkable effect from the interaction is 
that it gives rise to a negative thermal expansion coeffi- 
cient in some glasses, at low enough temperatures poj : 
The fluctuating entities in the lattice attract in the Van 
der Waals fashion, via exchanging phonons. This attrac- 
tion, counterbalanced by the materials stiffness, acts to 
partially contract the sample. The number of thermally 
excited transitions increases with the temperature and 
thereby enhances the degree of contraction. The effect 
is small, about 10~ 6 , but nevertheless observable. The 
often employed dimensionless parameter characterizing 
lattice non-linearity - the Gruneisen parameter - is usu- 
ally positive and of the order unity in crystals. It was 
found to be large and negative in many glasses at cryo- 
genic temperatures, see e.g. |38j . 

The direct phonon mediated TLS-TLS interaction goes 
as 1/r 3 , just like the usual dipole-dipole interaction, 
but is somewhat complicated by the tensorial form of 
the coupling (see e.g. [3(j for a discussion). If, how- 
ever, the transverse and longitudinal speeds of sound 
were equal, the interaction would be strictly dipole- 
dipole. It is therefore often convenient to assume the 
elastic interaction is indeed purely dipole-dipole result- 
ing in a small quantitative error. With this simplifi- 
cation, a "scalar" phononic Hamiltonian can be used: 
oKj,kl ( f ) ij ( t > kl \p c2 s{^4>) 2 i where is a scalar displace- 
ment field polarized in a single direction. The coupling 
will become gij4>ijO~ z — > (g\7(f>)o~ z . (The tensorial char- 
acter of the interaction may actually be important in 
the Gruneisen parameter context, see [2(j and below). 
The presence of an electric component to each transi- 
tion dipole moment clearly leads to another contribution 
to the total interaction. Since a detailed discussion of 
the interaction effects has already been given elsewhere 
[2fl ] , here we simply estimate the strength of the electric 
dipole-dipole coupling relative to the purely elastic coun- 
terpart and the rest follows in a straightforward fashion. 

The elastic dipole-dipole interaction is given by a sim- 
ple formula, see e.g. [20|. [T(ij ]: 

q 1 / a \ ^ 

Wdast - 9 3 ~ T 9 \ ~) > ( 21 ) 

where, for the sake of clarity, we eschew some numerical 
constants (these could be found in Ref. 20]) and have 
used Ea. (|17f> . The electric dipole-dipole interaction, by 
Ecis. lfT^ll and is, on the other hand: 

H clcct ^t 2 pcla 3 . (22) 

What is the relative value of the two interactions? A use- 
ful rule of thumb is that g ~ \jT g pc 2 s a? is of the order eV 



for all substances. In silica, for instance, pc 2 s a? is several 
eV, the Rydberg scale being a convenient (and physi- 
cally justified) landmark. The T g of silica is 1500 K, i.e. 
slightly larger than 0.1 eV. We therefore make an inter- 
esting observation that the electric dipole-dipole interac- 
tion can be comparable in magnitude to the elastic coun- 
terpart for polar enough substances. This is despite the 
relatively weak contribution (1% or less) of the coulom- 
bic forces to the apparent mechanical stiffness. We will 
speculate on the physical significance of this observation 
in the final Section of the article, while for now, we limit 
ourselves to a formal notion: The elastic dipole-dipole in- 
teraction is disadvantaged, compared to the electric coun- 
terpart, due to the large pc 2 s term in the denominator of 
Ea. (|2~T|l : Phonons are not true gauge particles. 

Leggett has emphasized that the dimensionless 
Gruneisen parameter varies "wildly" between different 
amorphous substances, in contrast with the nearly uni- 
versal / m fp/A ratio. Lubchenko and Wolynes have argued 
pfij . this stems from the Van der Waals attraction be- 
tween the tunneling centers, which is strongly enhanced 
by "Boson Peak" excitations. The total attractive inter- 
action consists of several contributions, is temperature 
dependent, and is expressed in terms of various combi- 
nations of the temperature, T g and the Debye tempera- 
ture. While it may be argued that there is an intrinsic 
upper bound on the value of £ (see below), there other- 
wise seems to be little intrinsic connection between the 
polar and structural characteristics of glasses, in general. 
As a result, the electrostatic interaction from Eq.JS^J) is 
expected also to contribute to the non-universality of the 
Gruneisen parameter. This notion is consistent with the 
sensitivity of the Gruneisen parameter in silica to the con- 
centration of polar impurities j^. Further, the magni- 
tude of the negative thermal expansivity is indeed larger 
in more polar mixtures according to Ref. [3S| . 

Symmetry vs. Transient Piezoelectricity. We have so 
far focused on the dipole moment of one of the structural 
states relative to the other, namely JV A/x^ = p^, since it 
is what determines the coupling of the transition to exter- 
nal field. We next inquire whether there is a correlation 
between the degree of polarity of a state with its abso- 
lute energy. For example, suppose for a moment that the 
lower energy state is completely non-polar, so that the 
two-level system dipole moment is completely due to the 
excited state. Such a situation could be exploited exper- 
imentally: One could supercool a liquid just enough so 
that it does not crystallize readily (the way glassblowers 
do), then expose the sample to strong electric field for 
a sufficient while, and then quench the sample below its 
T g . After that, remove the field. Clearly, the number 
of dipole moments along and opposite to the field direc- 
tion will differ. In other words, a number of dipoles will 
be actually lined up in a preferred direction, leading to 
a (weak) ferro-electric order. As a consequence of this, 
removal of the field in the procedure above should lead 
to sample's contraction (which is the sample's way of 
minimizing the ferroelectric energy.) Some residual po- 
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larization will appear as well. If the field is removed suffi- 
ciently fast, the sample will also heat up some. (This, in a 
sense, represents an antithesis to adiabatic demagnetiza- 
tion.) Such polarization will not take place if the degree 
of polarity is uncorrelated with the energy of a structural 
state, at least within the relevant energy range. 

While we can not, at present, rule out a priori the 
"pyroelectric" scenario above, it seems rather unlikely 
for the following reason. An amorphous sample, at least 
at long enough times, should have inversion symmetry, 
on average. With this symmetry present, no piezoelec- 
tricity, let alone pyroelectricity, can take place |24|. The 
question therefore is whether an external field breaks the 
inversion symmetry (locally) without inducing crystal- 
lization. Perhaps it will, on short enough time scales, 
before considerable aging takes place. An experimental 
study would settle this issue. At any rate, even if present, 
ferroelectric order would not affect the thermal expansion 
properties at low temperatures (c.f. the discussion of the 
Griineisen parameter in [20|h This attraction mecha- 
nism is temperature independent and simply contributes 
to the effective molecular field at each TLS site. 



IV. CLOSING REMARKS 

In this article, we have outlined the microscopic origin 
of the coupling of the intrinsic structural transitions in 
amorphous solids to electric fields. The coupling stems 
from rotation of the molecular bonds, within the region 
encompassing the transition, which generates a net elec- 
tric dipole. The molecular constituents of a glass, even 
if intrinsically non-polar, are strained due to disorder. 
Therefore small partial charges on each atom are ex- 
pected leading to the presence of electric dipole moments 
associated with individual bonds. 

A local structural transition occurs by moving a do- 
main wall through the region. The domain wall is a me- 
chanically strained region separating alternative struc- 
tural states. Such strained regions are frozen-in thermal 
fluctuations of the lattice. (Above T g , the life-time of 
such as frozen-in structure would not exceed the typical 
a-relaxation time in the liquid.) In an analogous fashion, 
the dipole moment can be thought of as frozen-in polar- 
ization fluctuations, or as the local piezoelectric response 
to local strains. 

The glass transition is driven by steric, i.e. mechanical 
interactions, not primarily by electrical ones, which is re- 
flected in the smallness of the partial charge, given in the 
theory by the dimensionless quantity £• So, for example, 
the elastic modulus can be represented as pc 2 s (l + ( 2 ), 
where C 2 S 0-01 is the contribution of the ionic forces to 
the overall material stiffness. (Here, we discriminate be- 
tween ionic (or hydrogen) bonds, as in NaCl or H2O, and 
covalent bonds, as in diamond.) In spite of its apparent 



small contribution to the material's structural integrity, 
the coulomb component could actually be comparable 
to the elastic component of the interaction between dis- 
tinct transitions. Note such interactions are unimportant 
as far as the identity of each transition is concerned, be- 
cause the corresponding regions are statistically far apart 
[2fj| . Nevertheless, this weak interaction may be viewed, 
loosely, as the successful attempt of the system to have 
avoided entropically costly local ferroelastic order. With 
this in mind, and the relatively strong coulomb interac- 
tion of structural fluctuations, the smallness of the partial 
charges in many glasses may simply reflect the fact that 
purely ionic compounds, such as sodium chloride or wa- 
ter, do not vitrify readily but instead, form low density 
lattices. (Note, ice does amorphize under high pressure, 
see e.g. |4J|.) Conversely, this imposes an upper bound 
on the £ 2 value. 

Consider now Ea. (|llfl . In order for the lattice to be 
stable, the A' matrix (with an ij pair considered a sin- 
gle index) should be positive definite. This means, if a 
particular A is negative, the corresponding (3 could be, 
in principle, quite large. This would imply large induced 
dipole moments. As just argued, the coulomb component 
is likely to be small in a glass, because the latter is rather 
homogeneous. Defects in crystals, however, can be more 
extended and can be highly anisotropic. In this regard we 
wish to mention the single chromophore studies of spec- 
tral drift in a Shpolskii matrix, by Bauer and Kador |4l| . 
These authors have seen a transition (presumably due 
to a structural defect) that generated a remarkably large 
dipole moment of 8 Debye. Note, the Shpolskii matrix is 
polycrystalline, not strictly amorphous. 

Finally we remark there is more to electrodynamics 
of amorphous solids than what could realistically be dis- 
cussed in this communication. Here we have analysed 
only the case of good electric insulators. Distinct, in- 
teresting phenomena take place in semiconductors and 
metallic glasses. As an example, let us mention an old 
experiment of Claytor and Sladek who found en- 
hanced ultrasonic attenuation in AS2S3 glass upon re- 
moval of electric field. This extra attenuation was greatly 
reduced by infrared radiation, which suggested that there 
is, we quote, "atomic relaxation accompanying electronic 
transition in gap states where injected carriers have been 
trapped" . We leave this for future work. 
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